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ABSTRACT 

We report on very high energy 7-observations with the MAGIC Telescope of the pulsar PSR 
B1951-I-32 and its associated nebula, CTB 80. Our data constrain the cutoff energy of the pulsar 
to be less than 32 GeV, assuming the pulsed 7-ray emission to be exponentially cut off. The upper 
limit on the flux of pulsed 7-ray emission above 75 GeV is 4.3 • 10~^^ photons cm~^ sec~^, and the 
upper limit on the flux of steady emission above 140 GeV is 1.5 • 10~^^ photons cm~^ sec~^. We discuss 
our results in the framework of recent model predictions and other studies. 

Subject headings: acceleration of particles — gamma rays: observations — - pulsars: individual (PSR 
B1951-t-32) — radiation mechanisms: nonthermal 
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1. INTRODUCTION 

It is currently believed that pulsars are among the few 
objects in our Galaxy that are candidate sources of ul- 
trarelativistic charged cosmic rays. Relativistic particles 
within the magnetosphere emit 7-rays at energies up to 
several GeV in various processes such as curvature radi- 
ation, synchrotron radiation and inverse Compton (IC) 
scattering. Thus, observations in the multi-GeV 7-ray 
domain allow one to study the acceleration sites in the 
magnetosphere of a pulsar. Predicted sites where particle 
acceleration can take place are, for example, above the 
polar cap of the neutron st ar (e.g. iHarding et akl 119781 : 
iDaughertv fc Harding! Il982f ) and in the so-called out er 
gap of the magnetosph ere (e.g. ICheng et al.l Il986al lbl: 
IChiang fc Romanil Il992f ). Furthermore, particle accel- 
eration can take place outside the magnetosphere in the 
region where the pulsar wind interacts with the interstel- 
lar medium. If electrons are accelerated in these shocks, 
they could give rise to IC-scattered photons from, for ex- 
ample, the cosmic microwav e background, synchrotron 
radiation, or thermal origi n (Ide Jager fc H arding" 19921: 
lAtovan fc Aharonianlll996l : iBediiarek fc Ba rtosik 2003.). 
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PSR B1951+32 was de tected first at radio frequencies 
bv lKulkarni et alJ ()1988l ). and is one of the six rotation- 
powered high energy pulsars whose GeV emissio n was de- 
tected by EGRET ( Ramana murthv et allll995f) . Among 
7-ray pulsars, PSR B1951+32 is the only source ob- 
served to emit up to 20 GeV with no cutoff being evi- 
dent in the differential energy spectrum. The spectrum 
shows a hard spectral index of 1.8 between 100 MeV 
and 20 GeV. The pulsar has an apparent high efficiency 
(~ 0.4%) of converting its rate of rotational energy loss, 
3.7 X 10"^® ergss^^, into 7-rays above 100 MeV (assum- 
ing a distance of 2.5 kpc to the pulsar). Moreover, the 
7-ray luminosity at ^ 10 GeV is comparable to that of 
the Crab pulsar (| Ramanamurthv et al.|[l995l ). 

AS inferred from its rotational parameters, the 
spin-down age of P SR 61951+32 is ~ lO^r 
(|Manchester et al-lbOOSD .^"^ that is, about 100 times older 
than the Cr ab pulsar. The m agnetic field strength of 
4.9 • 10" G (|Manchester et al.l 12005) is lower than that 
in most rotation-powered pulsars. Because of the lower 
magnetic field, curvature 7-rays emitted near the stel- 
lar surface, as predicted in polar-cap models, are less 
affected by magnetic pair production. Compared with 
younger, more strongly magnetized pulsars, the spectral 
cutoff energy is thereby shifted to higher energies, up 
to a fe w tens of Ge V (Harding (2001); Baring (2004); 
see also iBulik et al.l pOOOl ) for a discussion of low-field 
millisecond pulsars). 

On the contrary, if the 7-rays are emitted in the outer 
magnetosphere, as predicted in outer gap models, the 
potential drop in the outer gap of PSR B1951-I-32 is 
expected to be comparable to t h at of young pulsars 
(s ee eg. jUjoi [Zhang fc Cheng I (fl99l and eq. [2.1] 
of ICheng et all (ll986a()). Therefore, the cutoff energy, 
which reflects the maximum Lorentz factor of the elec- 
trons or positrons accelerate d in the outer gap, is ex- 
pected to be around 10 GeV ([Hirotanill2007h . Thus, fea- 
tures in the predicted spectral shape of weakly magne- 
tized pulsars at energies above 10 GeV are strongly de- 
pendent on the emission altitudes. In order to discrimi- 
nate between emission models, PSR B1951+32 is a prime 
candidate for observation by ground-based 7-ray detec- 
tors with low energy thresholds such as the imaging air 
Cherenkov telescope MAGIC. 

This pulsar is located in the core of the radio nebula 
CTB 80, which is thought to be physically associated 
with t he pulsar. In X-rays the nebula s hows a cometary 
shape (jMoon et al.ll2004l : iLi et all 120051 ). being confined 
by a bow shock that is produce d by the pulsar' s high 
proper motion (240 ± 4 km s ~^) (|Mighazzo et al.ir2002 ). 
iBednarek fc BartosikI (j2005b| ) predict an over-200GeV 
flux from the nebula at a level of ~ 4.4% of the Crab's 
flux, by assuming that high-energy leptons can accumu- 
late inside the well-localized nebula for long periods of 
time, as observed in the case of the Crab nebula. 

The current tightest constraint on the emission above 
100 GeV from the pulsa r and its nebula, obtain ed by the 
Whipple collaboration (jSrinivasan puts an 

upper limit of 75 GeV on the cutoff energy of the pulsed 
emission and an upper limit of 1.95 x 10~^^ cm~^s~^, on 
the steady emission above 260 Ge V. The latter is within 
a factor 2 of the prediction of IBednarek fc BartosikI 

See |http : / / www. atnf.csiro.au/ research /pulsar /psrcat | 



(|2005bH . 



TABLE 1 

Summary of the observations of PSR B1951- 
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Date 


Rate 


ON Time 


Extinctio: 


Q Extii 


iction Scatter 




(2006) 


(Hz) 


minutes 


(mag) 




(mag) 


Selected? 


Jul 4 


164 


130 


0.099 




0.017 




Jul 5 


164 


136 


0.100 




0.011 




Jul 6 


167 


105 


0.088 




0.014 


yes 


Jul 7 


176 


62 


0.091 




0.011 


yes 


Aug 3 


151 


95 


0.161 




0.009 


yea 


Aug 4 






0.266 




0.045 




Aug 23 


175 


168 


0.079 




0.017 


yea 


Aug 24 


158 


105 


0.088 




0.014 


yea 


Aug 25 


165 


138 


0.142 




0.029 


yea 


Aug 26 


135 


148 


0.168 




0.044 




Aug 27 


167 


124 


0.140 




0.042 


yea 


Aug 28 






0.249 




0.056 




Sep 13 


147 


83 








yes 


Sep 14 , , , 


139 


155 


0.105 




0.016 


yea 


Sep 15 


156 


102 


0.091 




0.017 


yes 


Sep 16 , , 


147 


125 


0.095 




0.013 


yes 


Sep 17 


149 


89 


0.094 




0.060 


yes 


Note. — Th 
the Carlsbcrg 
The cxtinctior 


Meridian 
I cocfficie 


on coeffici 
Telescope 


ents are take 
■, which is lo 
n effective w; 


:ivelength 


iblicly available 
the same site a 
of 625nm. 


data from 
s MAGIC. 



In this paper, we present upper limits on the cutoff 
energy of the pulsed emission from the pulsar, as well 
as on the steady and pulsed very high energy (VHE) 
fluxes from the region associated with the radio nebula, 
resulting from MAGIC telescope observations that were 
performed in 2006 July through September. The paper 
is structured as follows: After a short introduction to 
MAGIC and our data taking and analysis (fj2|), we re- 
port on our search for steady and pulsed emission from 
PSR B195H-32 (g]). We close with a discussion of the 
imphcations of our results (gj). 

2. OBSERVATIONS AND ANALYSIS 

The MAGIC (Major A tmospheric Ga mma Imaging 
Cherenkov) Telescope, sec 'Lorcn^ (j2004l ). is located on 
the Canary Island of La Palma (2200 m above sea 
level, 28.45°N, 17.54°W). MAGIC is currently the largest 
imaging atmospheric Cherenkov telescope (lACT), hav- 
ing a 17 m diameter tessellated reflector dish comprising 
964 0.5 X 0.5 diamond-milled aluminium mirrors. The 
faint Cherenkov light flashes produced by air showers are 
recorded by the telescope camera, which consists of 577 
photomultiplier tubes. Together with the current config- 
uration of the MAGIC camera w ith a trigger region of 
2.0° diameter (j Cortina et al.ll2005[ ). this results in a trig- 
ger collection area for 7-rays of about 10^ m^ at small 
zenith angles. The effective collection area depends on 
the analysis and is ~ 10"* m^ around 60 GeV and increases 
to > 6-10'' m^ beyond 200 GeV. At present, the minimum 
trigger energy is 50-60 GeV (at small zenith angles). The 
MAGIC Telescope is focused to 10 km distance — the 
most likely position for a 50 GeV air shower maximum. 
The accuracy in reconstructing the direction of incoming 
7-rays on an event by event basis (point spread func- 
tion), is about 0.1°, depending on energy and the chosen 
analysis method. A source with a 7-ray flux of ~ 2% 
that of the Crab Nebula and the same spectral slope can 
be detected by MAGIC above 200 GeV at a significance 
level of 5 (T within 50 hours. 
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Fig. 1. — Distribution of the parameter \alpha\ for events 
> 280 GeV. The distribution of OFF-source events (red) was nor- 
malized to the ON-source events (black) between 20° and 85°. 
An excess due to 7-rays from PSR B1951+32 is expected for 
alpha] < 7.5° {shaded region). 

PSR B1951+32 was observed with MAGIC for a total 
of 17 nights between 2006 July 4 and September 17. The 
observations were performed in the so-called ON/OFF 
mode; that is, PSR B1951-)-32 was observed by directly 
pointing to it (ON). Three nights were rejected because 
of unstable trigger rates due to bad weather. The back- 
ground was estimated by observing at the same range of 
zenith angle for 5.8 hr a suitable region in the sky where 
no 7-ray source is expected (OFF). In total, 30.7hr of 
data were processed. The zenith-angle range of the ob- 
servation was restricted to between 5° and 25°, guar- 
anteeing the lowest possible energy threshold. A sum- 
mary of the observations is given in Table [T] This table 
also includes the atmospheric extinction coefficients for 
all nights, provided by the Carlsberg Meridian Telescope, 
which is located at the same site as MAGIC. 

Following calibration of the data (|Gaug et al.l 120051 ) 
and a tail-cut image cleaning of the ev ents, a Hilla s 
parametrization algorithm was applied (|Hillasl I1985D . 
The tail cuts used in the image cleaning were 6 photoelec- 
trons for core pixels and 4 photoelectrons for boundary 
pixels. For the generation of sky maps, we used tail cuts 
of 10 and 5 photoelectrons. Additional suppression of 
pixels containing noise was achieved by requiring a nar- 
row time coincidence between adjacent pixels (~ 7nsec). 
The hadronic background was suppress ed with a mul- 
tivariate rnethod , the Random Forest (jBreima ^ [20011 : 
iBock et al.|[200l . which uses the Hillas parameters of an 
event to decide on its so-called hadronness. The power to 
suppress hadronic background is energy dependent and 
reduced for 7-ray energies below 150 GeV. As a conse- 
quence, the optimal cut in hadronness, which gives the 
highest rejection of background while retaining most of 
7-ray candidates, has to be independently determined 
for each energy region. For the analysis of the data 
presented here, we used an energy dependent hadron- 
ness cut, whose empirical parametrization was derived 
from Monte Carlo (MG) studies. An exception is the sky 
maps, for which a static hadronness cut was applied in 
the event selection. This is justified, as the maps were 
produced for energies above 200 GcV, where the depen- 
dence of the optimal hadronness cut on energy is small. 
The method of random forests is also used to estimate 



the energy of an event. Typically, energy resolutions of 
25% are achieved on an event-by-event basis. 

3. RESULTS 

3.1. Search for Steady Emission 

We searched for steady 7-ray emission of a point source 
from the direction of PSR B1951-I-32 with different anal- 
ysis thresholds between 140 GeV and 2.6 TeV. We define 
the analysis threshold as the peak of the energy distri- 
bution of MG events after cuts. Images of 7-rays from 
PSR B1951+32 point with their major axis to the cam- 
era center and thus appear as an excess at small values in 
the parameter "alpha". Alpha is the angle between the 
major axis of the shower image and the direction deter- 
mined by the image's center of gravity and the camera 
center. In Figure [U we show the distribution of |alpha| 
for events with energies > 280 GeV. An excess due to 7- 
ray emission from PSR B1951-I-32 should be visible in the 
figure for |alpha| < 7.5°. The results of this analysis and 
others with different analysis thresholds are summarized 
in Table [2] As no significant signal (> 5(t) from 7-rays 
was found, we calculated upper limits on the number of 
excess eve nts with a co n fidenc e level of 95% by using the 
method of iRolke et al.l ()2005( ). In the calculation of the 
limits, a systematic uncertainty on the flux of 30% was 
taken into account. The upper limits on excess events 
were converted into integral flux limits by assuming a 
spectral index of 2.6, which is similar to the spectral 
index of the predictions and other known pulsar wind 
nebulae (PWNs) such as the Crab Nebula. If a harder 
spectrum with index 2.0 is assumed, the flux limits in- 
crease by about 15% and they decrease by about 40% if 
a softer spectrum with index 4.0 is assumed. The inte- 
gral flux limits of 7-rays are shown in Figure [2] together 
with the meas urement of | Srinivasan et al.l ( )1997l ) and the 
predictions of lBednarek fc BartosikI (|2003l ). 

3.2. Search for "f- Ray Emission in the Vicinity of PSR 

B 195 1+32 

We explored the region in the sky around the position 
of the pulsar for a possible extended or displaced emission 
region of 7-rays. The latter is a likely scenario because of 
the high proper motion of the pulsar. For this stud y, we 
employed the DISP method of Fomin ct al.' ("1994 1) with 
a mo dified parametrization (Domingo-Santamari a et al.l 
l2005f ). which permits the reconstruction of the arrival 
direction of a > 100 GeV 7-ray with an accuracy of ~ 
0.1°. Sky maps were produced in different bins of energy. 
In none of the maps was 7-ray emission found within the 
reconstructed field of view, of ~ 0.6° radius. 

The map in Figure [3] {left) shows the significance cal- 
culated in bins of 0.1° x 0.1° for events with energies 
> 200 GeV. Figure [4] shows a map of the calculated up- 
per limits (95 % confidence level) on the integral flux for 
the same events. The acceptance of the MAGIC camera 
was modelled using the radial dependence of the back- 
ground rate in the camera after event selection. By com- 
paring with MC simulations, we confirmed for various 
angular distances from the camera center that the radial 
dependence of the background rate is compatible with 
the simulated 7-ray acceptance. 

Following our study, we can exclude steady 7-ray 
emission above 200 GeV at the level predicted by 
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TABLE 2 

Results of the Analysis Searching for Steady 7-Ray Emission from PSR B1951+32. 



Analysis Threshold 


Significance Upper Limit, 


Flux Upper Limit 


(GcV) 


ON Events OFF Events Excess Events {a) Excess Events (95% C.L.) 


(cm-^s-^) 



> 140 37869 37933 ±381 -64 -0.2 792 1.5 x 10"^^ 

> 280 3576 3740 ± 150 -164 -1.0 196 2.7 x 10"^^ 

> 530 712 777 ± 42 -65 -1.3 54 7.0 x 10"^^ 

> 800 232 231.5 ± 22 0.5 0.0 55 7.0 x 10""^^ 

> 1060 101 90.6 ± 14 10.4 0.6 45 5.8 x 10""^^ 

> 1400 58 49.5 ± 10.8 8.5 0.6 35 3.9 x 10"^^ 

> 2600 17 26 ± 10 -9 -0.9 14 2.5x10"^^ 
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Fig. 2. — Integral upper limits (95 % confidence level) on the 
steady 7-ray emission from the direction of PSR B1951+32. For 
comparison, the 7-ray flux of the Crab Nebula l lWagner et al.l2005l ) 
is also indicated. 



energy between 100 GeV and 2 TeV. To test for pe- 
riodicity, we ap plied the Pearson- test, the if- Test 
( de Jager et aLlil989) . and a test from lGregorv fc Loredol 
( 1992f ) (a Bayesian-Test). No signature of pulsed emis- 
sion was found in any of the energy intervals. As an ex- 
ample, we give the results of the ff-Test, which yielded 
significances of 0.3, 2.3, 0.6, 0.2 and 1.4 cr, respectively, 
with increasing energy. The corresponding 95% confi- 
dence level upper limits on pulsed emission are shown 
in Figure [5l T he limits were calculated from the results 
of the H-Test (|de Jageijll99l by assuming a duty cycle 
for the pulsed emission of 36%, which corresponds to the 
duty cycle of PSR B1951-I-3 2 at energies above 100 MeV 
(|Ramanamurthv et al.lll995l ). A spectral slope of 2.6 was 
assumed in the calculation of the upper flux limit. Note 
that these are upper limits in differenti al bins of energy, 
whereas the upper limits from Whipple (jSrinivasan et al] 
IT99I are integral ones, which were converted to differ- 
ential ones assuming a spectral shape of 2.6. 



iBednarek fc BartosikI (|2003f ). which we would have de- 
tected if (1) the emission were originating from within a 
circle of radius « 0.4° centered on the position of the pul- 
sar and (2) the apparent emission region was restricted 
to less than ^ 0.3° in diameter. 

3.3. Search for Pulsed Emission 

The time of each event (hereafter "arrival time" ) is de- 
rived from the time signal of a GPS-controUed rubidium 
clock with a precision of ~ 200 ns. Before we searched for 
pulsed emission from the pulsar, the arrival times were 
transformed to the barycenter of the solar system with 
the Tempo timing package by J. H. Taylor et al.^^ After- 
wards, the corrected arrival times tj were folded to the 
corresponding phase (pj of PSR B1951-I-32: 

(t>j = i^itj - to) + ]^v{tj - to f + ^Htj - to? 

where v and to are the values from a contemporary 
ephemeris provided by A. Lyne (2006, private commu- 
nication), which is listed in Table [3l The analysis chain 
that was set up to search for pulsed emission was previ- 
ously tested on data from an optical observation of the 
Crab pulsar with the central pixel of the MAGIC camera 
(|Lucarelli et al.|[2005[ ) . Details of the optical observation 
can be found in F. Lucarelli et al. (2007, in preparation). 

We performed a search for pulsed 7-ray emission from 
PSR B1951+32 in five differential bins of reconstructed 

■^^ See |http : / / www, atnf.csiro.au/research/pulsar /tempo | 



TABLE 3 
Ephemeris of PSR B1951+32 



Parameter Value 



Position epoch (JD) 2,450,228.4144 JD 

R.A 19^'52™58!27568995 

Decl 32°52'40.6824033" 

Pulsar Epoch (JD) 2,453,931.724208 JD 

V (Hz) 25.29516019929(63) 

V (Hz s-i) -3.72818(33) • IQ-^^ 

V (Hz s-2) -1.15(25) ■ 10-21 

Note. — From A. Lyne (2006, private communication). 
Uncertainties are given in parentheses. 

In a second analysis, we searched for pulsed emission 
by selecting events with a SIZE> 100 photoelectrons^^, 
that is, events with energies > 75GeV. Again, no hint of 
pulsed emission was found. The ff-Test yielded 1.4, and 
a test 7.2 with 11 degrees of freedom. The Bayesian 
test gave a probability for pulsed emission of 2.4 • lO"''. 

From the result of the ff-Test, we calculated an upper 
limit on the number of excess events (s. Table |4]), from 
which we derived an upper limit on the cutoff energy of 
the pulsed emission in the following way: The known 
spectrum of P SR B1951-t- 32 at GeV energies, measured 
by EGRET (,Fifi7rolfT995h . was multiplied by an expo- 

■^^ SIZE is the integrated intensity of a shower image after applied 
tail cuts in units of photoelectrons. It is also a good measure of 
the incident energy for shower impact parameters between ~ 50 to 
120 m. 
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Fig. 3. — Significance of VHE 7-ray emission from the region around PSR B1951+32. Left: Calculated significance of VHE 7-ray emission 
> 200 GeV in bins of 0.1° X 0.1°. O verlaid in black are contours from radio observations (Castelleti et al. 2003) and in white contours from 
IR observations (iFesen et al.lll98H) . Right: Distribution of significances. The distribution is compatible with that of randomly distributed 
data. 

nential cutoff and convolved with tlie effective collection 
area of the telescope. For a given cutoff energy, we then 
obtained the number of expected excess events by mul- 
tiplying the result with the dead-time-corrected obser- 
vation time. The upper limit on the cutoff energy was 
finally found by iteratively changing the cutoff energy 
until the number of expected excess events matched the 
upper limit on the number of pulsed excess events. With 
this procedure we obtained an upper limit on the cut- 
off energy of 32 GeV. The measured spectrum of PSR 
B1951-f 32 multiplied by an exponential cutoff of 32 GeV 
is shown in Figure^ (red curve). The analysis threshold, 
75 GeV, is marked by the red arrow in the figure. In case 
that the rollover of the 7-ray spectrum is superexponen- 
tial in shape, we constrain the cutoff energy to be below 
60 GeV. 

As a cross-check, the same analysis was repeated, 
this time by selecting all events with a SIZE < 
300 photoelectrons, that is, events with energies < 
180 GeV. The resulting pulse phase profile in Figure [S] 
shows no evidence for pulsed emission. From this anal- 
ysis, a slightly better upper limit on the cutoff energy 
of 28 GeV results. The analysis threshold, 60 GeV, was 
lower because events with a SIZE below 100 photoelec- 
trons were also included in the analysis. 

4. DISCUSSION 

Theoretical predictions and experimental evidence 
from lower energies had been quite favorable for a possi- 
ble detection of 7-ray emission from PSR B1951-I-32 or 
its nebula with MAGIC. Nevertheless, despite the higher 
sensitivity of this observation compared to previous ones, 
no 7-ray emission was detected. 

The upper limits in Figure [2] on the steady 7-ray emis- 
sion from the PWN surrounding PSR 1951+32 are below 
the 7-ray flux that was predicted by the time-dependent 
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Fig. 4. — Upper limits (95 % confidence level) on the integral 
7-ray emission above 200 GeV, calculated in bins of 0.05° X 0.05°. 



model of iBednarek fc BartosiS (|2003l l2005aD . Although 
their model takes into account the temporal evolution 
of the nebula (but not the spatial evolution), the ac- 
celeration of leptons and therefore also the equilibrium 
spectrum of leptons inside the nebula still depends on 
a few free parameters. These parameters, for example, 
the density of the medium surrounding the PWN, the 
acceleration efficiency of leptons, and the magnetization 
parameter of the pulsar wind at the shock region, are not 
well constrained by observations. 

Concerning the magnetization parameter, that is, the 
ratio of t he magn e tic en ergy flux to the particle en- 
ergy flux, iLi et al.l (|2005f) have recently estimated the 
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Fig. 5. — Results of the analysis in the search for pulsed emission 
from PSR B1951+32. Upper limits are given with a 95 % confi- 
dence level. The upper limit on the cutoff energy from Whipple is 
shown as the dot-dashed curve. The upper limit on the cutoff of 
32 GeV by MAGIC is shown as the solid red curve. The analysis 
threshold (75 GeV) is marked by the arrow on the horizontal axis. 



magnetic field strength of the compact X-ray nebula 
aromid PSR B1951+32 to b e ~ 300 fiG, which is l arger 
than the value assumed by iBednarek fc BartosiU . At 
the present time it is therefore clear that the value of 
the magnetization parameter a of the pulsar wind has 
to be much larger th an the value of a — 10^^, which 
IBednarek fc Bartosi k' assumed. As a result, the cooling 
of electrons by synchrotron radiation is faster and the IC 
7-ray flux is suppressed. Neverthel ess, a hadronic compo- 
nent, as predicted in so me models (|Bednarek fc BartosikI 
120031 : iHorns et al.1l2006f) . which would d ominate if the ac - 
celeration efficiency of leptons was low ('Bedn arekl[2007l ). 
would be below the sensitivity of our observation. 

Another aspect is that the model of 
IBednarek fc BartosikI deals with PWNe that are 
well confined by the external medium and pulsars 
that are, at most, moving slowly through the inter- 
stellar medium (the prototype of such a nebula is 
the Crab nebula). Only in such a scenario should a 
well-localized 7-ray source be expected, whereas when 
a pulsar is moving very fast, the 7-ray emission will be 
distributed over a larger volume. In the case of PSR 
B1951+32, which is movi ng with an apparent v elocity 
vpsR = 240 ± 40 km s^^ (iMidiazzo et al.l 1200^. the 
7-ray flux estimated by IBednarek fc BartosikI ( 20b5bD 
will be smeared over an area with a diameter d of at 
least 



WPSRTPSR ~ 5.3 • 10^^ cm w 0.5 



(1) 



assuming an age of the pulsar of Tpsn = 7 x 10'* yr and 
a distance of 2 kpc. Such an extended emission region 
reduces the detection probability with MAGIC. Apart 
from the pulsar's motion and the diffusion of leptons, 
their confinement and cooling, as well as their injec- 
tion rate into the interstellar medium over time, have 
to be taken into account. These parameters are un- 
known, and therefore, their influence on the extension 
of the 7-ray source is difficult to estimate. Assuming 
Bohm diffusion in a magnetic field of 3 x 10~^ G, one es- 
timates a diffusion length of 13 pc for 100 TeV leptons 
during the lifetime of the pulsar (jBednarek fc BartosikI 
l2005bi In this case the extension would marginally in- 
crease by 0.1° beyond what is expected from the mo- 



tion of the pulsar alone. If the magnetic field distri- 
bution is ordered, the diffusion can be faster and even 
anisotropic, leading to much larger emission regions. In 
this context it is interesting to note that extended TeV 
7-ray sources associated with displaced pulsars were re- 
cently detec ted by the H.E.S . S. Co llaboration (e.g., the 
Vela pulsar lAharonian et al.l (|2006f) and PSR B1823-13 
lAharonian et al.l (|2005[) V 




Fig. 6.— Pulse phase profile of PSR B1951+32 obtained after 
selecting events with SIZE < 300 photoelectrons. The shaded areas 
indicate the phase regions i n which PSR B 1951-f32 is emitting at 
GeV energies ( Ramanamur tliv et al.iri995l) . 

Considering the 7-ray emission from the pulsar, we 
constrain the cutoff of the pulsed emission to less than 
32 GeV if the cutoff is an exponential, which is appropri- 
ate when the 7-rays are emitted more than a few neutron 
star radii above the surface. If photons are emitted at 
lower altitudes, they are subject to magnetic pair pro- 
duction, resulting in a stronger (superexponential) at- 
tenuation of the energy spectrum. In the latter scenario 
we constrain the allowed range of cutoff energies to be 
< 60 GeV. Considering further that large uncertainties 
govern the last spectral point measured by EGRET, it 
follows that the allowed energy region where the cutoff 
resides can be constrained to lie somewhere between 10 
and 30 GeV (exponential cutoff) or up to 60 GeV (super- 
exponential cutoff). The narrow allowed range does not 
leave much freedom for models. This result and the up- 
per limits from the search in differential bins of energy 
are compared in Figure [5] with theoretical predictions 
from the polar-cap and the outer-gap model. In this fig- 
ure, the dotted line represents the polar-cap predictions 
from (Harding 2001), renormalized to the points of the 
EGRET spectrum. The thin solid l ine shows the spec- 
trum of the latest outer-gap model (|Hirotanill2007f l. 

In polar-cap models, the cutoff energy is determined 
by the attenuation of 7-rays due to magnetic pair pro- 
duction and hence by the emission altitude of 7-rays. 
As a consequence, the energy spectrum above the cut- 
off energy is superexponentially attenuated. If the emis- 
sion altitude in the polar cap model shown in Figure [5] 
changes from 1 to 2 stellar radii, the cutoff energy will 
increase from 20 GeV to 60 GeV, which is, according to 
our observations, the maximum allowed cutoff energy for 
a superexponentially shaped cutoff. On the contrary, in 
outer-gap models, the cutoff is determined by the maxi- 
mum Lorentz factor of the accelerated positrons and elec- 
trons. As a consequence, the cutoff of the 7-ray spectrum 
is smoother, resulting in an exponential cutoff. If the 
magnetic field lines near the light cylinder are straighter 
than assumed for the outer-gap spectrum in Figure |5l 
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2.4- 10-4 
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3.2 
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3388 


5.0 ■ 10-" 


10.7 


3.6 - 10-4 



the predicted flux below 60 GeV will increase. 

For more precise predictions of the cutoff energy in 
polar-cap models, multidimensional and self-consistent 
electrodynamics have to be examined from flrst princi- 
ples, whereas a three-dimensional magnetic field config- 
uration has to be investigated in the outer-gap model. 
Assuming that these improvements in theory will be 
achieved in the near future, measurements with higher 
statistics around 10 GeV, for example, by GLAST, or 
measurements by future ground-based experiments with 
lower thresholds than MAGIC, for example, MAGIC II 
or the Cherenkov Telescope Array, will be needed in or- 
der to distinguish between models. 

The predicted IC flux at TeV energies in the outer- 
gap model (Figure [5] solid black line) appears to be in- 
consistent with our upper limits. Nevertheless, it must 
be noted that the IC flux is obtained by assuming that 
all the magnetospheric soft photons illuminate the equa- 
torial region of the magnetosphere in which the gap- 
accelerated positrons are migrating outwards. Therefore, 
the predicted IC flux as a function of energy specifies 
an upper boundary to the possible pulsed TeV emis- 
sion. The open poloidal magnetic field lines could have a 
single-signed curvature within 1.8 light-cylinder radii, as 
the solution of the time-dependent force-free electrody- 



namics of an oblique rotator indicates (|Spitkovskvll2006() . 
If this is the case, soft photons emitted inside the light 
cylinder along the convex magnetic field lines will not 
efficiently illuminate the magnetic field lines, which are 
slightly convex even outside the light cylinder. As a re- 
sult, the predicted IC flux at TeV energies will be sig- 
nificantly reduced. This problem will be solved in future 
when the s elf-consistent gap electrodynamics ()Hirotanil 
I2OO6I l2007t ) and the three-dimensional force-free electro- 
dynamics are combined. 
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